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Aerothermodynamic Characteristics in the 
Hypersonic Continuum-Rarefied Transitional 

Regime 

Christopher E. Glass* and James X. Moss^ 

NASA Langley Research Center , Aerotherrriodynamics Branch , Hampton , K4 23681-2199 


Stagnation pressure and heating for cylinders and spheres representative of wing lead- 
ing edges and nose sections for small-scale space transportation vehicles are presented for 
various Knudsen numbers at Earth entry conditions that bridge from the free molecular 
to continuum regimes. CFD, DSMC, and collisionless DSMC are used to span Knud- 
sen numbers from 0.001 to 100 with the CFD/DSMC overlap at a Knudsen number of 
0.01 and the DSMC/collisionless DSMC overlap at 10. The present results include den- 
sity contours, density and temperature along the stagnation streamline for all cases, and 
stagnation pressure and heat transfer coefficient as a function of Knudsen number. 


Nomenclature 

Ch Heat transfer coefficient, q' f / 

C p Pressure coefficient, ( p - Poo)/^Po cA ^ 

D Diameter, m 

K rip Knudsen number based on diameter, D 

p Pressure. Pa 

q n Heat transfer rate, W/nr 

R Radius, m 

T Temperature, K 

V Velocity, m/s 

x, y, z Cartesian coordinate system, m 
Z Geometric altitude, km 

p Density, kg/m 3 

Subscripts 

av Overall, associated with all modes 

Iran Translational 

rot Rotational 

vib Vibrational 

ic Wall 

oc Free stream 

Introduction 

A CCURATE and credible prediction of three- 
dimensional leading edge aerothermodynamic 
characteristics in the transitional flow regime at and 
near rarefied flight conditions bridging from contin- 
uum to rarefied flow is a challenge. For small-scale 
space transportation vehicles, including small exper- 
imental X- vehicles, crew transfer vehicles (CTY), or 
crew and cargo transfer vehicles (CCTV), and future 
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air breathing vehicles having extremely sharp load- 
ing edges, rarefaction phenomena air extremely im- 
portant as maximum heating may occur on leading 
edges in the rarefied, as opposed to continuum, regime. 
Also, heating loads may be different than inferred 
from continuum approaches. These mispredictions in 
the aerothermal environment may have corresponding 
negative effects in the I PS material selection and siz- 
ing. 

For predicting the flight aerothermodynamies in the 
transitional flow regime, numerical techniques are re- 
lied on to a greater extent to define leading edge flow 
field and surface effects than experimental techniques. 
Few experimental rarefied flow facilities are avail- 
able 1 to test models representative of complex, three- 
dimensional space transportation vehicles. Also, for 
most experimental wind tunnel model scales of these 
vehicles, leading edge dimensions do not lend to accu- 
rate discrete measurements because the sensor size is of 
the same magnitude as the leading edge and, for a sin- 
gle data point, integrates the loads over a significant 
portion of the leading edge. Three-dimensional heat 
conduction affects small model leading edge measure- 
ments because it diffuses the flow induced heating and 
reduces the 1 maximum heating value. Also, global lead- 
ing edge measurement techniques that visualize the 
surface, such as pressure sensitive paint, infrared ther- 
mography, or thermographic phosphors, suffer from 
rapid pixel- to- pixel viewing angle changes, insufficient 
pixel coverage, and three-dimensional conduction as 
with discrete heat transfer measurements. 

Numerical techniques also have their limitations, 
which must be addressed. Commonly used numerical 
techniques to provide 1 aerothermodynamic estimates in 
the transitional flow regime at and near rarefied flow 
conditions are computational fluid dynamics 1 (CFD) 
and direct simulation Monte Carlo 2 (DSMC). In gen- 
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oral. CFD is applicable to the continuum regime, al- 
though with some approximations to the near surface 
How. it can provide solutions into the transitional flow 
regime. Additionally. CFD offers a complication of 
finely gridding the flow about leading edges, which pro- 
duce large flow gradients, and joining these regions to 
coarser gridding over acreage regions of the body. 

Of available numerical techniques, DSMC offers the 
most accurate and credible solution method for the 
rarefied and near rarefied condition. However, be- 
cause the method is probabilistic. i.e., it tracks move- 
ments and collisions of simulated molecules and sta- 
tistically samples the positions and surface interac- 
tions. it requires long run times on a computer. In 
addition, as the flow condition approaches the con- 
tinuum from the rarefied regime, the DSMC method 
requires many more simulated molecules and the at- 
tendant larger computer memory and longer CPU run 
times. Yet this computer hardware restriction has re- 
cently been diminished with the introduction of larger 
multi-processor computers and DSMC software writ- 
ten specifically for three-dimensional problems using 
parallel processing algorithms to take advantage of 
multi-processor computer architectures. 

Motivation 

The intent of this study is to employ a three- 
dimensional DSMC algorithm' 3 to assess its usefulness 
in providing three-dimensional leading edge flow field 
and surface solutions for flow at and near rarefied 
flight conditions with application to current and fu- 
ture aerospace vehicles. The particular issues that 
are investigated in this study are verifying the three- 
dimensional DSMC results by comparison with a vali- 
dated 2-D/axisymmetric DSMC algorithm 1 and CFD 
algorithm. 0 determining the bridging relationship of 
the stagnation point values over a range of An p from 
continuum to free molecular, and providing guidance 
for the transition from CFD to DSMC and from DSMC 
to collision less DSMC. 

By verifying the three-dimensional results in the 
transitional flow regime, greater acceptance will be 
gained in applying the DSMC technique to obtain 
aerothermodynamic definition on leading edges, where 
the technique more closely represents the flow physics 
than the the current CFD options. As computer ca- 
pability increases both in speed and memory, three- 
dimensional DSMC will become 1 the defacto method 
for near rarefied transitional flow regime solutions, and 
verification of the technique, such as performed in this 
study, will prove useful in providing confidence in ex- 
tending the applicable Knp range of DSMC. 

In addition, results of this study will be useful in 
providing applicability limits of CFD. DSMC, and col- 
lisionless DSMC techniques. In the present study, 
CFD is performed for Knp = 0.001 and 0.01, DSMC 
for 0.01 < Knp < 10. and collisionless DSMC for 



Fig. 1 Similarity between simple geometric shapes 
and wing leading edges and nose sections of small- 
scale space transportation vehicles. 

Knp = 10 and 100. At Knp — 0.01, CFD and 
DSMC results overlap, and at Knp ~ 10. DSMC 
and collisionless DSMC.' results overlap. The flow 
regime changes from continuum to transitional and 
from transitional to free molecular, respectively, near 
these Knudsen numbers, and by performing overlap 
calculations at these conditions and analyzing the re- 
sults. the suitability of each technique at the overlap 
can be assessed. 

Methodology 

For this study, the strategy chosen is to apply 2- 
D and 3-D CFD and DSMC techniques to simple 
cylindrical and spherical geometries at various flight 
conditions that span 0.001 < Knp < 100. Shown in 
Figure 1 are tw r o flight vehicles that are designed to 
travel through this transitional regime. Also shown on 
the figure is the geometric similarity between a cylin- 
der and sphere, which are presently studied, and wing 
leading edges and nose sections. 

A diameter of D = 0.1 m w r as chosen for this 
study because it is close 1 to the wing diameter of the 1 
small-scale space transportation vehicles shown in Fig- 
ure 1. Further, the w'all is assumed noncatalytic with 
a T w — 300 K to produce large temperature variation 
through the shock layer and to exercise the techniques 
to a greater extent near the wall. For all cases, a con- 
stant velocity of — G000 m/s is specified, the 

air properties and composition w r as obtained from the 
1976 U.S. Standard Atmosphere. 6 and reacting gas 
chemistry is utilized. Table 1 shows free stream con- 
ditions at various altitudes taken from the standard 
atmosphere. 

A sketch of the coordinate system for the cylinder 
and sphere numerical solutions is presented in Fig- 
ure 2. Note that the cylinder and sphere are centered 
in the x, y . z system with the free stream flow' proceed- 
ing in the x-direction. Also, in the results presented 
subsequently, the coordinates are normalized by the 
radius, R. 
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Table 1 Free stream conditions for the A rip range. 


Kn D 

Z. km 

Too, K 

Pro- kg/m J 

P oc ^ 

0.001 

51.9 

269.2 

8.130 x 10- ' 

62.83 

0.01 

70.5 

219.2 

7.722 > 

; lO " 5 

4.837 

0.1 

84.0 

188.6 

8.051 > 

c 10 6 

0.4360 

1 

98.0 

191.7 

8.071 > 

c 10 ~ 7 

0.04506 

10 

111.4 

256.3 

7.657 > 

c 10~ 8 

0.006009 

100 

132.0 

488.7 

6.904 > 

c lO" 9 

0.001109 



Fig. 2 Cylinder and sphere coordinate system. 

As mentioned previously, the molecular analysis 
is performed by the DSMC Analysis Code (DAC), 3 
which is a three-dimensional code with both scalar 
and parallel processing options. DAC has achieved 
a significant acceptance for rarefied three-dimensional 
applications. The parallel version is currently em- 
ployed to produce DSMC solutions over the majority 
of the Krin range. Collisionless DSMC provides so- 
lutions in the high Knp region: G2 code' 1 results are 
utilized to anchor and verify DAC results throughout 
the the Knp range; and GASP 5 CFD provides contin- 
uum solutions and bridges to the near continuum DAC 
results. Table 2 presents the numerical technique and 


Table 2 Numerical technique and code for cylinder 
and sphere results. 


Knp 

CFD 

DSMC 

Collision less DSMC 

0.001 

GASP 



0.01 

GASP 

DAC.G2 


0.1 


DAC 


1.0 


DAC.G2 


10 


DAC.G2 

DAC 

100 



DAC 

specific 

computer 

code used to provide 1 solutions for 


the various Knudsen numbers. As shown in the table, 
the number of flow solutions is 11 each for the cylinder 


and sphere. 

The DAC code requires a triangulated surface grid 
of the body of interest, which is then positioned into 
a Cartesian flow field grid, for the DSMC solution. In 
this study, the 2-D. axisymmetric, and 3-D options of 
the DAC code' wen 1 exercised. Therefore, tlirn? sur- 
face grids were built for the DSMC portion of the 
study, one for the cylinder and two for the sphere. In 
addition, the DSMC solutions included the wake re- 
gion so that the entire surface geometry was modeled. 
The surface. 1 grid for the 2-D cylinder consisted of 360 
triangles to completely enclose the geometry, and the 
sphere was composed of 12,363 surface triangles. Note 
that a 5° slice of the sphere with 1 46 triangular ele- 
ments is the geometry for the axisymmetric DAC case 
(Knp = 0.01). 

Both the cylinder and sphere CFD solutions wore 
only about the windward portion of the geometry and 
the outflow boundary condition used is an extrapola- 
tion from the interior, which is appropriate because the 
flow exiting the grid is supersonic. Because the solu- 
tion was 2-D, tin 1 grid used for the cylinder consists of 
257x361x2 points with 257 points away from the wall. 
For the sphere, the grid is composed of two blocks. The 
block with the stagnation region is a 37x37x81 point 
block with 81 points away from the wall. Tin 1 stagna- 
tion region block was constructed with no singularity 
at the stagnation point so that difficulties associated 
with solutions about a singular point were avoided. 
The second block of the sphere flow field patched the 
stagnation region with the outflow boundary. For both 
geometries, the grid normal to the wall is closely clus- 
tered to resolve the large gradients at the stagnation 
point and spacing proceeds away from the wall with 
an exponential factor of 3 = 1.001. 

Fully reacting chemistry models are used for both 
the DSMC and CFD solutions. DAC solutions employ 
translational, rotational, vibrational, and chemical en- 
ergy transfer modes between six molecules and atoms 
(N 2 , 0 2j NO, N, (). and Ar). The present applica- 
tion of DAC uses rotational and vibrational relaxation 
numbers of 5 and 50, respectively. GASP CFD solu- 
tions utilize the non-equilibrium, finite-rate, 5-specie 
(N 2 . 0 2 , NO, N, and 0) chemistry model of Kang and 
Dunn with vibrational energy taken into account for 
the three diatomic molecules. 

Results and Discussion 

This section presents results from the 22 flow so- 
lutions that are given in Table 2. Discussed in the 
first section is numerical uncertainty of the molecular 
and continuum techniques. Next, results consisting of 
flow field density contours, temperature and density 
variation along the stagnation streamline from tilt' free 
stream condition to the wall, maximum atomic oxygen 
and nitrogen along the streamline, temperature (over- 
all, translational, rotational, and vibrational) variation 
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along the stagnation streamline for the K n p — 0.01 
case, and the stagnation pressure and heating are 
given. The contours will show the effect of rarefac- 
tion of the overall flow field. Presented next, tem- 
perature and density along the stagnation streamline 
more clearly show the effect of varying Knudsen num- 
ber on the shock layer structure. Also, along the 
stagnation streamline, the maximum mole fraction of 
atomic oxygen and nitrogen and the overall, transla- 
tional, rotational, and vibrational temperature mode 
distribution will allow an assessment of the degree of 
chemical and thermal nonequilibrium of the flow. Fi- 
nally. the stagnation pressure and heating as a function 
Knudsen number will exhibit the degree of agreement 
between the various techniques and the bridging of 
these surface properties from the continuum to the rar- 
efied regime. 

Properties presented in this section arc normalized 
by an appropriate factor. Temperature and density are 
normalized by the free stream values given in Table 1. 
The normalization factors for C p and Cm are given in 
Tabled. 


Table 3 Normalization factors for C p and Ch. 


Kn n 

P>D 


1/2 p 


1 /2poo l 

0.001 

6.283 x 1U 1 

1.463 

x 10' 

8.780 x 10 7 

0.01 

4.837 x 

: 10° 

1.390 x 10 3 

8.340 x 10 6 

0.1 

4.360 x 

io-' 

1.449 

x 10- 

8.693 x 10 5 

1 

4.506 x 

10~ 2 

1.453 

x to 1 

8.717 x 10' 

10 

6.009 x 

10“ 3 

1.378 

x 10° 

8.269 x 10 3 

100 

1.109 x 

10 :i 

1.243 : 

t 10 1 

7.436 x 10- 


Numerical Uncertainty 

A major source of numerical uncertainty associated 
with applying both the DSMC and CFD techniques 
is the user defined grid; therefore, a brief discussion 
of this source of uncertainty specific to the present 
study is presented. As mentioned previously, the DAG 
DSMC requires a triangulated surface geometry de- 
scription, which is placed in a Cartesian grid and the 
simulation is performed within the grid. One of the 
primary uncertainty sources is the grid in a high gra- 
dient region such as a shock wave or flow near the 
body. The DAC series reduces this uncertainty source 
by including a preprocessing code that produces an 
adapted Cartesian grid. Application of the regridding 
process is controlled bv the user. For this study, mul- 
tiple adaptations were performed to ensure that the 
grid in regions of large gradients and near the surface 
were dimensioned to less than the mean free path res- 
olution. Sub-mean free path gridding is required to 
resolve completely the surface quantities with DSMC 
(For example, see Moss, et al. 7 ). 

For CFD. flow field convergence, which is a function 
of both the solution L l > norm and the grid density, 
affects the solution uncertainty. The CFD solutions 



Fig. 3 Cylinder flow field density contours for 
I\nn = 0.001 from CFD. 


of the present study were advanced through a coarse, 
medium, and fine grid to accelerate convergence and 
allow a comparison of the solution at the various grid 
densities. After the solution was converged on the 
coarse grid, it was interpolated onto the finer grid 
and the solution converged again. The grid density 
change from the coarse to medium grid is accomplished 
by increasing the number of computational cells by a 
factor of 2 all directions; likewise, the change from 
a medium to fine grid is produced by increasing the 
number of cells by a factor of 2 in all computational 
directions. For the present CFD cases, when the L 2 
norm decreased between 5 to 7 orders of magnitude 
and the L 2 norm, surface pressure, and heating became 
steady, the solution was assumed converged, in addi- 
tion, when converged solutions on the medium grid, 
were compared with the fine grid solution, the stag- 
nation pressure difference was less than about 2% and 
the stagnation heating less than about 5%. 

Density Contours 

Density contours for the Ktip 0.001, 0.01. 10, 
and 100 cases are presented in this section for both 
the cylinder and sphere. These Knudsen number cases 
were chosen because they represent the end points of 
study {Ktip 0.001 and 100) and the crossover points 
where two techniques were applied (Ktip 0.01 and 
10 ). 

Cylinder 

Shown as Figures 3. 4, 5. and 6 are the flow field con- 
tours for a cylinder at the Ktip - 0.001. 0.01. 10, and 
100, respectively, and the free stream conditions given 
in Table 1. For the continuum flow case (Figure 3), 
the CFD solution of the flow field shows a sharp bow 
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Fig. 4 Cylinder flow field density contours for 
Knp = 0.01 from DSMC and CFD. 


Fig. 5 Cylinder flow field density contours for 
K rip = 10 from DSMC and collisionless DSMC. 

shock with a high density stagnation region. The con- 
tours also show that as the flow proceeds both above 
and below the cylinder, the density decreases because, 
away from the stagnation region, the flow is allowed 
to expand. 

Shown as Figure 4 are the density contours for the 
Knp - 0.01 condition. At this Knudsen number, both 
DSMC and CFD solutions were obtained because the 
condition is near the transition from continuum to 
transitional flow. The flow field shown on the upper 
portion of the figure is taken from the DSMC solution 
and the one on the lower portion is taken from the CFD 

5 


Fig. 6 Cylinder flow field density contours for 
Knp = 100 from collisionless DSMC. 

solution. These density contours are very similar with 
the exception of the diffuse nature of the DSMC bow 
shock leading edge, which is most apparent where the 
two contours meet along the y/R = 0 line. Because 
the DSMC technique more closely simulates the prob- 
abilistic nature of the flow physics at this condition, 
the upper contour probably better describes the flow. 
Also, note that the DSMC solution models the flow 
about the cylinder including the wake region shadowed 
by the cylinder. A comparison of the DSMC and CFD 
flow field at the CFD outflow boundary (x/R = 0) 
shows little difference, which lends credence to not 
modeling the wake region with the CFD. 

The next contours, shown in Figure 5, are at the 
Knp - 10 flow condition, which is near the transi- 
tion between the transitional and free molecular flow 
regimes. The comparison shown on the figure is be- 
tween a DSMC and collisionless DSMC solution. Be- 
cause free molecular flow is collision less, the DSMC 
code with the collision routine disabled properly rep- 
resents a free molecular condition. 

Figure 5 shows that for either case, no distinct bow 
shock is present, only a diffuse region with a grad^ 
ual density increase as the flow approaches the cylin- 
der. The collisionless DSMC solution for this Knudsen 
number has a higher density region further forward of 
the cylinder than the DSMC solution. The 2-D nature 
of the cylinder flow with no relief in the /-direction 
docs not allow molecules to travel non-preferentially 
after a collision. Moreover, this effect shows that for 
the 2-D Knp - 10 flow condition, collisionless DSMC 1 
is not applicable to define this flow. 

The final cylinder flow contour shown is Figure 6 
for Knp = 100. This solution was produced using 
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Fig. 7 Sphere flow field density contours for 
Knn = 0.001 from CFD. 



Fig. 8 Sphere flow field density contours for 
Knn = 0.01 from DSMC and CFD. 

collisionless DSMC because the mean free path. A, of 
the free stream is about 10 m. and a free molecular 
assumption for the 0.1 m diameter cylinder at this 
condition is appropriate. No prominent bow shock is 
shown in the figure, only a diffuse high density region 
extending forward of the cylinder. 

Sphere 

Similar to the cylinder contours. Figures 7, 8, 
9. and 10 show density contours for the sphere at 
Knp = 0.001. 0.01. 10, and 100, respectively. The 
free stream conditions for these cases are taken from 



Fig. 9 Sphere flow field density contours for 
Knn = 10 from DSMC and collisionless DSMC. 



Fig. 10 Sphere flow field density contours for 
Knn = 100 from collisionless DSMC. 


Table 1 and the flow field density shown in the figures 
is normalized by the free stream density given in the 
table. 

Figure 7 shows the CFD flow field for the most con- 
tinuum condition (Knn = 0.001). Note that a sharp 
bow shock is produced in front of the sphere as the 
flow encounters it. When compared with the cylinder 
flow field at the same Knudsen number (Figure 3), the 
sphere bow wave is closer to body because the spheri- 
cal geometry provides a 3-D flow relief, thus allowing 
the closer bow shock. 

For the Knp — 0.01 case. CFD and DSMC were 
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employed to produce the How Held solutions shown 
in Figure 8. The flow Held at this free stream con- 
dition is near the transition between continuum and 
transitional flow. Note that a I) AC DSMC axisym- 
metric solution was required for this case to provide 
the grid density needed to accurately describe the near 
wall flow gradients with sub mean free path resolution. 
As with the cylinder results, DSMC density contours 
are on the upper half of the figure and CFD contours 
are on the' lower half. Tlu* DSMC bow shock is more 
diffuse at its leading edge than the CFD because, as 
discussed previously, the DSMC provides a better rep- 
resentation of this flow regime than the CFD. Again, 
the' DSMC flow field includes the wake portion behind 
the.' spin to and. as was seen with the* cylinder flow field, 
the CFD contours at the outflow boundary are similar 
to the DSMC contours at the some location. 

The' next density contours for the sphere are shown 
in Figure 9 for Knp 10. which is near the' Knudsen 
number between the transitional and free molecular 
flow regimes. Like the cylinder contours shown in Fig- 
ure 5, the contours on the upper half of the figure are 
from the' DSMC solution and the 1 collisionless DSMC 
contours are on the lower half of the figure. The bow 
shock is not distinguished with a sharp gradient asso- 
ciated with continuum flow, and the density contours 
in front of the sphere gradually rise as the.' flow moves 
close t to the sphere. Unlike the 1 cylinder results where 
the flow field is disturbed further upstream, a compar- 
ison of the two sphere density contours on Figure 9 
shows similarity. The 3-D flow relief of the sphere 
simulation allows the collisionless technique to better 
model this condition than the 2-D collisionless sim- 
ulation used for the cylinder. Therefore, the use of 
a collisionless simulation at this Knudsen number is 
better justified for a 3-D simulation. 

The Knp - 100 density contour figure for the 

sphere (Figure 10) is the final one presented in this 
section. Collisionless DSMC is employed to produce a 
solution for this free molecular flow field. Compared 
to the 2-D cylinder density contours (See Figure G.), 
the sphere 1 contours do not extend as far forward as 
with the cylinder. As discussed previously, the nature 
of the sphere simulation allows the molecules, which 
encounter the sphere, to reflect into a 3-D space, thus 
lessing the forward reflected molecules and better mod- 
eling the flow. 

Stagnation Streamline Results 

The variation of normalized temperature and den- 
sity along the stagnation streamline in the flow direc- 
tion for each of the GASP CFD, DAC DSMC, and 
collisionless DAC DSMC cases, given in Table 2, is 
presented in this section. The cylinder and sphere 
results are grouped from the continuum to the most 
rarefied condition in the following order: first, temper- 
ature and density for the CFD cases (Knp ~ 0.001 



Fig. 11 Cylinder stagnation streamline flow field 
properties for I\nn — 0.001 and I\np = 0.01 from 
CFD. 

and 0.01): second, temperature and density for the 
CFD and DSMC cases near the transition between 
the continuum and transitional regimes (Knp = 0.01); 
third, temperature and density for the DSMC cases 
(Knp =- 0.01. 0.1. 1. and 10); next, temperature 
and density for DSMC and collisionless DSMC near 
the transition between transitional and free molecular 
flow (Knp = 10); then, temperature and density for 
collisionless DSMC (Knp = 10 and 100); the maxi- 
mum atomic oxygen and nitrogen: and finally, overall, 
translational, rotational, and vibrational temperature 
distributions for Knp = 0.1. 

For the cylinder CFD cases, Figure 11 shows, for 
increasing Knudsen number, or decreasing free stream 
density, shock standoff distance increases and the tem- 
perature and density gradient at the leading edge of 
the bow shock decreases. Also, for the Knp = 0.01 
case, the initial temperature rise before encountering 
the bow shock occurs earlier, or further from the body, 
than for the Knp ~ 0.001 case 1 showing the ('fleet of 
a less continuum like flow. In addition, the maximum 
normalized temperature for Knp = 0.01 case is larger 
than for Knp — 0.001. 

The sphere CFD results at the same conditions are 
given as Figure 1 12. Although the sphere flow field 
properties show the same trends as those 1 of the cylin- 
der, they occur much closer to the surface because the 
3-D relief of the flow about the sphere. Also, note that 
for both the cylinder and sphere cases, a sharp gradi- 
ent in both temperature and density occurs near the 
wall because of the imposed 300 K wall temperature. 

Next, at Knp = 0.01, a comparison between the 
CFD and DSMC stagnation streamline temperature 
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x/R 


Fig. 12 Sphere stagnation streamline flow field 
properties for Kn n — 0.001 and Knn — 0.01 from 
CFD. 



Fig. 13 Cylinder stagnation streamline flow field 
properties for Knn — 0.01 from CFD and DSMC. 

and density for the cylinder, shown in Figure 13, and 
for the sphere, shown in Figure 14, is presented. The 
figures show that temperature and density gradients 
for DSMC are less than those for CFD. which was also 
shown in the density contour comparisons at the same 
Knudsen number (See, for example, Figures 4 and 8 
for the cylinder and sphere, respectively.). Note, how- 
ever. that the comparison for the cylinder is better 
than that for the sphere the maximum temperature 
rise and the temperature and density wall gradient are 
closer. Even though the CFD and DSMC cylinder re- 



x/R 


Fig. 14 Sphere stagnation streamline flow field 
properties for Knn = 0.01 from CFD and DSMC. 



Fig. 15 Cylinder stagnation streamline flow field 
properties for Knn = 0.01 to Knn = 10 from 
DSMC. 

suits are in better agreement than those for the sphere, 
in general, the proper treatment of the solution at this 
condition is bv performing DSMC, which allows for a 
more* realistic shock wave structure to evolve* as com- 
pared to the CFD calculation. 

Figures 15 and 16 show the normalized tempera- 
ture and density along the stagnation streamline for 
the cylinder and sphere, respectively, over the Knud- 
sen number range 0.01 < Kn p <10 as predicted by 
the DSMC technique. For both geometries, as the 
flow is initially affected by the body, the tempera* 
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Fig. 16 Sphere stagnation streamline flow field 
properties for I\nn = 0.01 to Knp — 10 from 
DSMC. 



Fig. 17 Cylinder stagnation streamline flow field 
properties for Kilo — 10 from DSMC and collision- 
less DSMC. 

ture and density increase, but the gradient decrease's 
as the Knudsen number increases. Also, the maximum 
overall temperature for Knp =0.1, shown in both fig- 
ures, is greater than for the other Knudsen numbers. 
Therefore, the thermal nonequilibrium aspects of the 
Knp = 0.1 case will be examined in greater detail 
subsequently. 

A comparison between DSMC and collisionless 
DSMC for a cylinder at Knp = 10 is shown next as 
Figure 17. As discussed previously with the compari- 
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Fig. 18 Sphere stagnation streamline flow field 
properties for Knp = 10 from DSMC and colli- 
sionless DSMC. 

son of Figure 5. the collisionless DSMC allows reflected 
molecules to return upstream without colliding with 
downstream moving molecules, which causes higher 
upstream density contours than the DSMC case. A 
comparison of the temperature distributions shows the 
collisionless DSMC temperature is much greater than 
that of the DSMC (See Figure 17.). This is because the 
collisionless overall temperature shown on the figure is 
only composed of the translational component, a gross 
measure of the average molecular speed, which for 
this case, has contributions from molecules streaming 
in both directions; therefore, showing the collisionless 
assumption for the cylinder at this condition is not 
appropriate. 

The comparison for the sphere at the same condi- 
tions is shown as Figure 18. Lnlike the cylinder, the 
sphere streamline temperature results for the two tech- 
niques are in much better agreement. In addition, the 
sphere density distributions are in good agreement, 
whereas, the cylinder density distributions shown in 
Figure 17 arc not similar. Overall, the cylinder and 
sphere results show the difference between imposing a 
2-D collisionless DSMC flow assumption, which is very 
restrictive, and 3-D collisionless DSMC at this Knud- 
sen number. 

The final comparison along the stagnation line is 
between the collisionless DSMC results at Knp = 10 
and 100, which are shown on Figures 19 and 20, 
respectively, for the cylinder and sphere. Although, for 
the cylindrical flow field, the Knp - 10 collisionless 
DSMC" has been shown not to be appropriate', it is in- 
cluded as a reference to compare with the Knp = 100 
case. The cylinder temperature and density distri- 
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Fig. 19 Cylinder stagnation streamline flow field 
properties for Knp = 10 and Knp = 100 from col- 
lisionless DSMC. 



Fig. 20 Sphere stagnation streamline flow field 
properties for Knp = 10 and Knn = 100 from col- 
lisionless DSMC. 



Fig. 21 Maximum atomic oxygen and nitrogen 
mole fraction along stagnation streamline. 

bution comparison shows that both Knudsen number 
cases exhibit the same behavior over the x/R range 
shown; that is, the overall temperature drops as the 
flow encounters the cold wall; Also, the density distri- 
butions are nearly the same. 

The sphere results, shown in Figure 20, however, 
have a different temperature distribution, which rises 
to a maximum value close to the sphere, then drops 
as the flow encounters the cold wall. This rise and fall 
is characteristic of three dimensional scattering of the 
molecules from the sphere surface. 

The maximum mole fraction of atomic oxygen and 
nitrogen along the stagnation streamline is presented 
next as Figure 21. This figure is included to show the 
degree of chemical nonequilibrium along the stagna- 
tion streamline for Knudsen numbers over the range 
0.001 < Knp < 1. The figure includes both the 
CFD and DSMC results for the cylinder and sphere. 
Line segments are also included between the calculated 
points in the figure only to distinguish tin 1 cylinder 
from the sphere results, not to provide an interpolation 
between the data points. The trend shown in the figure 
is that a maximum 0 and N in the shock layer occurs 
for Knp = 0.001, with a significant 0 mole fraction 
concentration (greater than 30%). Note that the O 
mole fraction is greater than the N because 0 > dissoci- 
ates at a lower energy level than N 2 . Also, the cylinder 
flow field shows a greater degree of dissociation than 
the sphere because the cylinder shock layer is thicker, 
which provides a greater residence time for the flow’ 
dissociation to proceed. In addition, the figure shows 
that for Knp >0.1. chemical dissociation becomes 
less significant in the flow* calculation. Therefore, it 
is essential to include a chemistry model in a flight 
environment stagnation flow calculation for the more 
continuum-like portion of the transitional regime, that 
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Fig. 22 Cylinder stagnation streamline flow field 
temperature for Knp — 0.1 from DSMC. 

is, for the range Knp < 1. 

As previously indicated, the maximum temperature 
along the stagnation streamline for the Knp = 0.1 
case, shown for the cylinder and sphere in Figures 15 
and 1G, respectively, was greater than for the other 
Knudsen number cases and may show the effects of 
thermal nonequilibrium; therefore, the overall, transla- 
tional, rotational, and vibrational temperature compo- 
nents for these two cases are presented and discussed. 
Also note, as shown in Figure 21 and discussed previ- 
ously, chemical nonequilibrium effects are present, but 
not significant for this case. 

Results for the cylinder are shown in Figure 22. and 
for the sphere, in Figure 23. The extent of thermal 
nonequilibrium is indicated by the temperature pro- 
files for each energy mode. Note that for this case, as 
shown in the figures, all three energy modes are ac- 
tivated inside the shock layer as the flow approaches 
the body. The energy distribution of the flow is re- 
lated to the translational and internal (rotational and 
vibrational) temperatures. The cylinder flow has a 
greater rotational and vibrational temperature than 
the sphere flow (Compare Figures 22 and 23.), per- 
haps because the thicker cylinder shock layer allows a 
greater degree of equilibration than the sphere shock 
layer, which may, in part, be chemically frozen and 
have a higher degree of thermal nonequilibrium be- 
cause of the thinner sphere shock layer. In any event, 
to accurately predict this flow field, a technique, which 
accounts for local thermal nonequilibrium effects of the 
flow, is needed. In addition, Cuda and Moss 8 and 
Moss, et al. 9 found significant thermal nonequilibrium 
effects from bodies at flight conditions near the same 
altitude as for the present Knp — 0.1 case. 
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Fig. 23 Sphere stagnation streamline flow field 
temperature for I\up — 0.1 from DSMC. 

Although not presented graphically, an examination 
of the various temperature modes for the cylinder and 
sphere at the Knp — 0.01. 1, and 10 conditions was 
performed. At Knp = 0.01, the flow in the shock layer 
was nearly in thermal equilibrium, and at the higher 
altitudes {Knp 1 and 10). the internal energy modes 
were essentially frozen at the free stream condition and 
were not in thermal equilibrium. 

Stagnation Point Results 

Presented are the stagnation point values from the 
22 cases that are shown on Table 2. These single point 
values of pressure and heating are given in coefficient 
form as C p and Cp with the normalization constants 
taken from Table 3. The Knudsen number range' of 
the data is 0.001 < Knp < 100, which spans from 
the continuum to the free molecular regime for the 
cylinder and sphere conditions, and defines the C p and 
Ch bridging relations for the two geometries over the 
Knudsen number range. Figure 24 shows the cylinder 
results, and Figure 25 shows the sphere results. 

Some of the results from the 3-D DSMC DAC code 
of LeBeau 3 that are presented in these figures an' an- 
chored by results from the 2-D/axisvmmetric G2 code 
of Bird. 2 A comparison of pressure and heat trans- 
fer coefficient is provided at Knp — 0.01, 1, and 
10 for the cylinder and sphere, respectively, in Fig- 
ure's 24 and 25. The pressure coefficient comparisons 
at Knp = 0.01 are in excellent agreement and at the 
other two Knudsen numbers are in good agreement for 
both the cylinder and sphere. The heat transfer coef- 
ficient. however, is in excellent agreement for all of the 
cases. Based on this, the 3-D DAC DSMC is shown to 
provide comparable stagnation pressure and heating 
results as the 2-D/axisymmetric G2 code. 
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Cw at various Knn from different techniques. 
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Ch at various Knn from different techniques. 


The pressure coefficient for both the cylinder and 
sphere, shown in Figures 24 and 25. respectively, 
starts at a continuum value of about C p = 1.9 at 
Klip 0.001 and increases to a value of about 2.1 at 
the free molecular condition (Knn — 100). There is a 
noticeable difference between the cylinder and sphere 
C p at the end points of the data where the cylinder re- 
sults are closer to C p — 2 than the sphere results. Also, 
there is a dip in the cylinder C p value at Knp = 0.01. 
Moss, et al 10 show a dip in drag coefficient at about 
the same Knudsen number for the Pathfinder aeroshcll 
shape. In addition, the collisionless DSMC cylinder re- 
sult at I\np 10 is higher than the DSMC result, but 
for the sphere, they arc. 1 nearly the same. As discussed 
above and shown by comparisons in Figures 17 and 18, 
the reason the sphere DSMC and collisionless DSMC 


results more closely match is that the 3-D flow field of 
the sphere is better modeled than that of the cylinder. 

Heat transfer coefficient for the cylinder and sphere 
is now discussed. For both geometries, the range 
of heat transfer coefficient is 0.1 < O/ £ 1 from 
the continuum (Knp 0.001) to the free molecular 
(Knp = 100) regime, respectively. In addition, the 
cylinder Ch rise lags the sphere Ch rise with increas- 
ing Knp by about 1/2 an order of magnitude. This 
may be because* tlu* sphere shock layer is thinner and 
closer to the* body with attendant higher gradients. 
These higher gradients in flow properties promote rar- 
efaction effects more quickly than with the thicker 
cylinder shock layer with longer flow resident time 
scab' for energy relaxation. At Knp - 0.01, the* point 
where both CFD to DSMC solutions were obtained, 
the CFD Ch is larger for both the sphere and cylin- 
der; however, the sphere CFD result is greater than 
the cylinder result, which is close to the DSMC result. 
This may show that CFD can provide an adequate 
cylinder stagnation heating estimate at Knp = 0.01, 
but cannot for the sphere. In addition, previous stud- 
ies that provide CFD and DSMC solutions as a func- 
tion of Knudsen number show that the CFD predic- 
tions are greater than DSMC predictions at the contin- 
uum transition (For example, see Cuda and Moss 8 and 
Gupta. 11 ). Based on the results at Knp = 0.01 where 
both CFD to DSMC techniques were implemented, 
CFD may bo used to obtain stagnation values for the 
cylinder, but, for the sphere, DSMC is needed. 

At Knp - 10, which is the 1 condition that both 
the DSMC and collisionless DSMC solutions were ob- 
tained, the coliisionless DSMC cylinder and sphere 
results are greater than the DSMC with the differ- 
ence between the cylinder results greater than with 
sphere results, which have little difference. Because 
both sphere stagnation heating results are close at the 
Knp =10 condition and the cylinder results are not, 
a change in solution methods from DSMC to collision- 
less DSMC is better justified at this condition for the 
sphere than for the cylinder. Also, for a 3-D blunt 
configuration, Moss, et al. 10 show that free molecular 
flow begins at a Knudsen number of order 10. 

Concluding Remarks 

CF D, DSMC, and collisionless DSMC were utilized 
to provide solutions for a cylindrical and spherical ge- 
ometry at reentry flight conditions, which span the 
Knudsen number range 0.001 < Knp < 100. In par- 
ticular, state-of-the-art computer programs an 1 exer- 
cised for this problem to assess the applicability. Den- 
sity contour comparisons, stagnation point streamline 
temperature and density distributions, and stagnation 
point values are shown for most of the cases. In ad- 
dition, a comparison between 2-D/axisymmetric G2 
results and 3-D DAC results is provided to anchor the 
3-D DAC results; The results at the stagnation point 
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compare well and show the 3-D DAC DSMC is equiv- 
alent to the G2 DSMC. 

At Knp = 0.01, near the continuum change to the 
transitional regime, CFD and DSMC solutions were 
obtained. For both the cylinder and sphere flow. CFD 
and DSMC flow fields differ, with the DSMC shock 
layer extending further from the surface. However, un- 
like the sphere results, the CFD and DSMC cylinder 
stagnation point values are nearly the same. Similarly, 
at Kni) = 10, the DSMC and collisionless DSMC are 
compared near the change from the transitional to the 
free molecular regimes. Neither the flow field nor the 
surfa.ee quantities for the cylinder compare well be- 
tween the DSMC and collisionless DSMC. However, 
the spherical flow results, which properly simulate the 
3-D nature of the Knp = 10 condition, match more 
closely. 

Chemical and thermal nonequilibriuin effects were* 
identified by examining the atomic oxygen and nitro- 
gen concentration and the translational, rotational, 
and vibrational temperature components along the 
stagnation streamline for various eases. For the flight 
condition in the transitional regime, chemical nonequi- 
librium effects are significant at Knp = 0.001 and 
should be included for calculations where hup < 1. 
Thermal nonequilibriuin was also shown to be impor- 
tant, especially for h up > 0.1. 

In addition, the bridging of the cylindrical and 
spherical stagnation point pressure and heating over 
the 0.001 < Knp < 100 is presented. For the con- 
ditions of the study, pressure coefficient, varies little 
from a value* of C p = 2; however, heat transfer coeffi- 
cient varies from about 0.1 at the continuum limit to 
1 at the* free molecular limit. 
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